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REVIEW OF BLOOD VESSELS VELOCITY ESTIMATIONS 


IRINA — ANDRA TACHE ! 


Rezumat. Bolile vasculare se reflecta prin tulburari ale curgelui sdngelui cat si 
modificari morfologice. Imagistica medicala digitala a_ revolutionat elaborarea 
diagnosticul si tratamentul acestor boli. Extragerea maximului de informatie despre 
fluxul sanguin din achizitiile medicale imagistice are o mare importanta pentru medicii 
din zilele noastre. Lucrarea prezinta o trecere in revista a celor mai semnificative metode 
de extragere a vitezei sdngelui din imaginile medicale. Angiografia cu raze X este 
standardul in diagnosticarea bolilor vasculare si poate fi utilizata cu succes impreuna cu 
metodele de estimare a vitezei sangelui prezentate in aceasta lucrare. 


Abstract. The vascular diseases are reflected by impairments in both morphology and 
hemodynamics. Digital medical imaging has revolutionized the diagnosis and treatment 
of these diseases. The extraction of the blood flow information from the daily clinical 
medical imaging acquisitions is one of the greatest importance for the physicians 
nowadays. The paper presents a review of the most noticeable methods of extracting 
blood velocity from medical images. X-ray angiography is the gold standard of the 
vascular diseases’ diagnosis, and it can be used successfully along with blood velocity 
estimation methods presented in this paper. 
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1. Introduction 


The cardiovascular diseases are the leading cause of death worldwide, followed 
by cancer as stated by World Health Organization. Medical imaging is intensively 
used in the diagnostics of vascular diseases in the past tens of years. The gold 
standard for the diagnosis of the coronary artery diseases remains the X-ray 
angiography. 

The main goal is the extraction of the maximum knowledge from a medical 
images. There are efforts in extracting the velocity from fluoroscopic angiography 
which offers important information about the propagation of the contrast agent 
into blood stream. 


2. Blood Velocity Methods 
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2.1 Bolus tracking methods 


The transit time method is used to compute the time density curve between two 
regions of interest (ROIs) along the same vessel, proximal and distal, and estimate 
the transport time. 

The approximated distance between the proximal and distal locations can be 
determined as Euclidian distance between the centerline pixels. Based on the 
assumption that the contrast agent perfectly mixes with the blood and it travels 
with the same speed into the vessels, the mean blood velocity can be estimated as 
the ratio between vessel length and transit time between the ROIs. This method 
was successfully applied for non-pulsatile flows of the cerebral and coronary 
vessels [1]. 


2.2 Velocity determination from distance density curves 


The method is based on obtaining, shifting, and matching the distance density 
curves at certain moments of time. Seifalian et al. in [2] defined the spatial shift 
Lo as the best match of the curves. Finally, the fluid velocity is equaled with the 
ratio between the best shift do and time interval between the two profiles: 


Lo 
= — 1 
v(t) re (1) 


Hawkes et al. in [3] concluded that if the contrast agent bolus moves completely 
from the tagged vessel area between two consecutive frames, the method would 
fail. 

To reduce the noise sensitivity, Rhode et al. in [4] developed a blood flow 
waveform shape model using the principal component analysis for fitting the 
concentration distance curves. An improvement of this method is found in [5]. 
After this step the optimal shift between successive fitted distance density curves 
is computed. They concluded that the concentration-distance curve matching 
algorithm is the most successful from all the distance density curve techniques. 


2.3 Velocity determination from flow map 


Colchester et al in [6] built the iso-concentration flow map or the parametric 
image to overcome the limitations of the bolus tracking methods. The parametric 
image contains information about the variations of the vessel densities in time 
along its length. Therefore, its columns contain the series of density distance 
curves and its rows the time density curves. The velocity is computed as the slope 
of the iso-concentration contours in the parametric image [7]: 


ey 8 (2) 
v=a, 7 tan 
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Where L is the distance along the vessel, t - time and 0@ - inclination of the iso- 
contours. 

Brunt et al in [8] has applied this method to Digital Subtraction Angiography 
(DSA) for determining the carotid blood velocity and Bateman et al in [9] 
combined them with the parametric imaging for quantifying the blood flow rates 
from the volume flow and total cardiac output. 


2.4 Indicator-dilution method 


Cunningham et al. in [10] adapted the Stewart—Hamilton dye dilution method for 
stenosis assessment by computing the axial velocity - v(x) as: 


v(x) = (3) 


A(x) fy c(x,t)dt 
Where A(x) is the surface area of the stenosis vessel. 


2.5 Inverted continuity equation method 


This method uses the continuity equation of the local mass conservation law 

which is written for moving fluids as: 

© +V-v) =F-V2D (4) 

Where D is the contrast agent density, v is the blood velocity, F is the diffusion 

coefficient. 

Assuming a negligible diffusion, the coefficient F is set to zero, and the above 

equation is written for the single dimensional case [11]: 
ons ) + v(x,t)- a v =f (5) 

Where the velocity v(x,t) is determined as the ratio between the partial 

differentials of density in function of time and distance: 

dD/ at 

~ @D/ ax (©) 

An additional prerequisite must be taken for the term dD/ 0x. It was observed 

[12] that for small values it can introduce noises. Therefore, it must be high 

enough in order not to increase the velocity errors. Also, the applicability of this 

method for overall velocity measures is not recommendable [1]. Nevertheless, the 

assumption of zero-diffusion condition is valid only for the contrast agent which 

travels near the tip of the bolus leading edge. 


v(x,t) = 
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2.6 Optical flow methods 


They are intensively used in computer vision and it is dedicated to compute the 
velocity of rigid objects in every pixel from temporal image series. 
Imbert et al. in [13] developed an optical flow method for symmetric axis and 
parabolic flow patterns: 

aD x—x,)*\ aD 

Vmax "3 (1 : z)+5 
x I t 

Where Vmax is the maximum axial velocity for parabolic flows, x- is the axial 
position of the vessel’s centerline and r is the vessel radius with symmetric axis. It 
was successfully applied for human femoral artery quantification. 
A weighted optical flow algorithm is proposed by Rhode et al. in [4] by averaging 
the estimated velocity along a vessel with the magnitude weightings of the spatial 
derivative. Another approach of the optical flow algorithm was adapted to DSA 
images by Bonnefous et al in [14]. 


=0 (7) 


2.7 Fluid continuity method 


The one-dimensional equation of the mass conservation law to volume of fluids 
computed from two axial locations xy,and xy, which delimitate the vessel 
segment can be expressed as in [1]: 


D(xy,t) D(xy,, 
[ “[D(x, ty,) — Deere + af es ) — eat = 0 (8) 


Where A is the cross-sectional area and Q is the flow rate. This mathematical 
expression considers a negligible diffusion and a uniform contrast mixing. The 
flow rates are fitted using least squares algorithm for overlapping or variable sizes 
of vessels’ segments and computation of the global flow rate Q between the times 
ty, and ty, 


2.8 Model based approach 


The computational flow dynamics was a recently adapted to medical imaging. It 
models sufficiently accurate the blood vessels flow, for normal or abnormal 
geometries. This makes it desirable for clinical practice, but its drawback is the 
high computational power. A solution could be the parallelization on graphical 
parallel units or cloud computing. 

Wachter in [15] proposed an overall model for predicting the iodine concentration 
after the injection into a vessel. The main hemodynamic components taken into 
consideration are the shape of the waveform model, the contrast agent injection 
model, the mixing of blood and contrast agent model at the injection site and the 
contrast agent propagation model. 
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2.9 Velocity estimation after 3D reconstruction 


This is a research topic which appeared in scientific papers shortly after the 
rotational angiography [16]. In these scientific papers a rotational C-arm 
angiograph is used for generating a more accurate 3D reconstruction of the vessels 
and the 2D fluoroscopic angiography projections for extracting the vessels’ 
contrast agent irrigation patterns. Different methods for mapping the 3D and 2D 
information [17] are proposed for these image sequences. 

The spatial velocity is determined based on the bolus transportation time 
estimated from fluoroscopic angiography and the absolute length and radius of the 
vessel segment of the re-projected centerline obtained after 3D reconstruction of 
the rotational angiography [16]. This information can be further used as the 
boundary conditions in the computational fluid dynamics models [17]. 

The rotational transit time [15] is estimated directly from the 3D image 
acquisition to extract the flow parameters from the flow map. Both are further 
used as parameters in the identification of the physical model of the blood flow 
and the contrast agent transportation [18]. 


Conclusions 


The physiological and structural information of the vascular system are evaluated 
by the physicians in dedicated clinics using medical imaging for diagnosis and 
even treatment of the coronary artery diseases. 

The paper presents a brief review of the velocity estimation methods varying from 
analyzing the propagation of the contrast agent into the blood stream, until 
modeling the blood flow dynamics. Most of these methods can be applied for 
medical imaging techniques, such as X-ray monoplane or rotational angiography, 
magnetic resonance angiography, etc. 


Acknowledgment 


This work was supported by the grant for project Processing of CT medical 
images for the identification and evaluation of lung damage due to COVID19 
using fractal analysis and artificial intelligence techniques, Scientific session of 
young researchers in the competition AOSR-TEAMS, edition 2022-2023. 


REFERENCES 


[1] Shpilfoygel, S. D., Close, R. A., Valentino, D. J., Duckwiler, G. R., X-ray videodensi- 
tometric methods for blood flow and velocity measurement: A critical review of litera-ture, Med. 


Phys., vol. 27, pp. 2008-2023 (2000) 


10 Irina - Andra Tache 


[2] Seifalian, A. M. Hawkes, D. J., Colchester, A. C. F., Hobbs K. E. F, A new algorithm for 
deriving pulsatile blood flow waveforms tested using simulated dynamic angio-graphic data, 
Journal of Neuroradiology, Volume 31, Issue 3, pp 263-269 (1989) 

[3] Hawkes, D. J., Seifalian, A. M., Colchester, A. C., Iqbal, N., Hardingham, C. R., Bladin, C. 
F. and Hobbs, K. E. Validation of volume blood flow measurements using three-dimensional 
distance-concentration functions derived from digital x-ray angio-grams, Investigative Radiology 
Volume 29(4), pp. 434-442 (1994) 

[4] Rhode, K, Ennew, G., Lambrou, T., Seifalian, A.,. Hawkes D., In-Vitro Validation of a 
Novel Model-Based Approach to the Measurement of Arterial Blood Flow Waveforms from 
Dynamic Digital X-ray Images, Medical Image Computing and Computer-Assisted Intervention — 
MICCAT 2001, Lecture Notes in Computer Science Volume 2208, pp 291-300 (2001) 

[5] Rhode, K. S., Lambrou, T., Hawkes, D. J., Seifalian, A. M. Novel approaches to the 
measurement of arterial blood flow from dynamic digital X-ray images, IEEE Transac-tions on 
Medical Imaging 24(4), pp. 500-513. (2005) 

[6] Colchester, A., Brunt, J. Measurement of vessel calibre and volume blood flow by dy-namic 
quantitative digital angiography, Journal of Cerebral Blood Flow & Metabolism 3, pp. 640-641 
(1983) 

[7] Brunt, J. N., Wicks, D. A., Hawkes, D. J., Seifalian, A. M., du Boulay, G. H., Colches-ter, 
A. F. and Wallis, A. The measurement of blood flow waveforms from X-ray angi-ography. Part |: 
Principles of the method and preliminary validation, Proceedings of the Institution of Mechanical 
Engineers 206(2), pp. 73-85 (1992) 

[8] Brunt, J. N., Wicks, D. A., Hawkes, D. J., Seifalian, A. M., du Boulay, G. H., Colches-ter, 
A. F. and Wallis, A. The measurement of blood flow waveforms from X-ray angi-ography. Part 2: 
Application to video recordings of digital subtraction angiography, Proceedings of the Institution 
of Mechanical Engineers 206(2), pp. 87-91 (1992) 

[9] Bateman W.A., Kruger, R.A., Blood flow measurement using digital angiography and 
parametric imaging, Medical Physics Volume 11, pp. 153 — 157 (1984) 

[10] Cunningham I., S. Yamada, B. Hobbs, A. Fenster, Arterial flow characterization with a 
photodiode array-based imaging system,’’ Med. Phys. 16, pp. 179-187 (1989) 

[11] Efron, U., Price, R. R., Smith, C. W. and Brill, A. B. Method to determine the instanta- 
neous blood-flow using cine- or video-densitometric data, Proceedings of the SPIE 143, pp. 154 - 
161 (1978) 

[12] Horm K., Schunk, B., Determining optical flow, Artif. Intel. 17, pp. 185—203 


Review of Blood Vessels Velocity Estimations 11 


[13] Imbert B., J. Meunier, R. Mongrain, J. Hudon, M. Bertrand, (1997), Stenosis parameter 
assessment from contrast medium tracking in cineangiography with an optical flow method, in 
Medical Imaging 1997, Physiology and Function from Multidimensional Images, edited by A. 
Hoffman, Proc. SPIE 3034, 631-640 (1981) 

[14] Bonnefous O., Pereira, V.M., Ouared, R., Brina O., Aerts, H., Hermans, R., van Nijnat-ten, 
F., Stawiaski, J., Ruijters D., Quantification of arterial flow using digital subtraction angiography, 
Medical Physics, Volume 39 (10), pp. 6264-6275 (2012) 

[15] Wachter I., 3D Reconstruction of Cerebral Blood Flow and Vessel Morphology from X- 
Ray Rotational Angiography, Phd thesis, Department of Medical Physics and Bioen-gineering, 
University College London (2009) 

[16] Castro M. A., Putman C. M., and Cebral J. R., “Patient-Specific Computational Model-ing 
of Cerebral Aneurysms with Multiple Avenues of Flow From 3D Rotational Angi-ography 
Images,” Academic Radiology, vol. 13, no. 7, pp. 811-821, (2006) 

[17] Schmitt H, Grass M, Rasche V, Schramm O, Haehnel S, Sartor K., An X-ray-based method 
for the determination of the contrast agent propagation in 3-D vessel structures, IEEE Transactions 
Medical Imaging, Volume 21(3), pp. 251-62. (2002) 

[18] Waechter-Stehle, I., Groth, A., Bruijns, T., Brinac, O., Ruefenacht, D. A., Kulcsar, Z., 
Mendes-Pereirac, V., Perrenc, F., Hawkese, D. J., Weese, J., Model-based blood flow 
quantification from DSA: quantitative evaluation on patient data and comparison with TCCD, 


Proc. of SPIE Medical Imaging, Volume 8314, pp. 83144R.1-17 (2012) 


